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We present a method of calculating the electric charge density of glass and silica surfaces in contact 
with aqueous electrolytes for two cases of practical relevance that are not amenable to standard 
techniques: surfaces of low specific area at low ionic strength and surfaces interacting strongly with 
a second anionic surface. 



I. INTRODUCTION 

Ionization processes in aqueous solutions have long 
been a point of central interest to physical chemistry. 
Much progress has been made in recent years in under- 
standing the charging properties of such different entities 
as small molecules, polyelectrolytes and various kinds of 
interfaces pj. 

In this context, silica and silicate glass surfaces im- 
mersed in water are known to acquire a negative surface 
charge density, primarily through the dissociation of ter- 
minal silanol groups. The degree of dissociation and thus 
the surface charge density results from an equilibrium 
between countcrions at the glass surface and free ions in 
the bulk electrolyte. Experimentally, this type of equilib- 
rium and its dependence on the solution conditions can 
be studied by potentiometric acid-base titrations on col- 
loidal dispersions of non-porous silica particles Q . This 
technique actually measures the volume concentration of 
protons transferred between the surfaces and the solu- 
tion. In order for the surfaces to accommodate a suffi- 
cient amount of charge, the electrostatic interaction be- 
tween the surface sites must be screened at least partially 
by added salt ions, and/or the available surface area must 
be large. 

These constraints can be relaxed to some degree by re- 
sorting to alternative techniques like microelectrophore- 
sis streaming potential measurements JB , conduc- 
tomctry || and electroacoustic methods Mg]. All of 
these methods, however, rely heavily on approximate 
models for electrostatic or hydrodynamic processes in the 
interfacial region, introducing uncertainties that are diffi- 
cult to estimate. We are not aware of any way to measure 
directly the surface charge of silica in a solution of very 
low ionic strength. 

Theoretical studies of low ionic strength solutions typ- 
ically deal with dense colloidal and macroionic systems 
and consider the regime of "no salt" , "low salt" , or "coun- 
terions only" . Here the defining assumption is that the 
overall ionic strength is due predominantly to the par- 
ticles or macroions and the compensating counterions in 
solution, whereas the concentration of any additional ions 
is negligible. For colloidal dispersions this assumption is 
again legitimate if the specific surface area carrying the 
colloidal charge is large. 



A series of recent experiments has spurred interest in 
the charge on glass and silica surfaces of low specific area 
in pure water, i.e. systems for which the usual picture 
of the "no salt" regime does not apply. For example, 
interaction measurements using digital video microscopy 
and optical trapping suggest that highly charged latex 
spheres may experience an anomalous long-ranged at- 
traction when confined by charged glass walls [p[-p^[, 
contrary to the predictions of Poisson-Boltzmann theory 
||l3| [l5|| . In one particular case the attraction ap- 
pears to result from a hydrodynamic interaction driven 
by the spheres' electrostatic repulsion from a nearby wall 
fl6f . This explanation hinges on the heretofore untested 
assumption that the glass wall carries an effective charge 
density of —2000 ± 200 e/fim 2 , where e is the elementary 
charge. Such hydrodynamic coupling cannot explain the 
like-charge attractions measured for spheres confined be- 
tween two charged glass walls [|9|— pTTf] . How the walls' 
charge influence colloidal electrostatic interactions is not 
yet resolved, in part because of open questions regard- 
ing the charging state of the glass. We recently reported 
that the pair interaction of silica spheres remains mono- 
tonically repulsive even in the presence of a charged glass 
wall J!?]]. The spheres' effective surface charge density of 
— 700 ± 150 e//j,m 2 extracted from these measurements is 
considerably smaller than the value posited in Ref. |T^] 
for a compact glass surface. 

Because a silica surface's charge density depends on 
the local chemical environment, it necessarily varies with 
proximity to other charge-carrying surfaces. The in- 
terpretation of typical particle deposition experiments 
|I8| and force measurements by atomic force microscopy 
(AFM) [|l9|j2C| l or total internal reflection microscopy 
(TIRM) p,l\ for example, is complicated by the fact that 
the charge densities of the substrate and the probe are 
a function of their separation, a phenomenon known as 
"charge regulation" . Since local properties of the en- 
closed solution rather than bulk properties determine the 
charging state, a naive use of charging data from bulk 
measurements can lead to errors. 

In the following section we discuss how the experimen- 
tally supported I-pK Basic Stern Model for silica sur- 
faces may be used to calculate the elusive charge of glass 
plates and strongly diluted silica particles in deionized 
water. In the remainder of this paper we take advantage 
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of a recently proposed theoretical treatment of charge 
regulation (2^j2j| to discuss the the charge of a silica-like 
surface in close proximity to a second anionic surface, 
which will be chosen, in view of the most common appli- 
cations, as either of the same type or of constant charge 
(like sulfate latex) or of carboxylic nature (like carboxyl 
latex and many biological surfaces). 



II. EFFECTIVE CHARGE OF GLASS AND 
SILICA IN DEIONIZED SOLUTIONS 

The principal mechanism by which glass and silica sur- 
faces acquire a charge in contact with water is the disso- 
ciation of silanol groups H 



22| , giving rise to the surface charge density 
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Further protonation of the uncharged group is expected 
only under extremely acidic conditions [ p4|]25fl and will 
be disregarded. Similarly, we will not take into account 
the protonation of doubly coordinated Si2 — O groups as 
these are generally considered inert [ p4| ]. 

In addition to the hydronium or other counterions dis- 
sociated from the surface, the bulk electrolyte also in- 
cludes ions due to the autodissociation of water; the lat- 
ter can dominate the concentration of mobile ions if the 
ratio of surface area to solution volume is exceedingly 
small. Under these conditions, the charging state of the 
surface is controlled by the ionic strength and pH of the 
bulk electrolyte, just as in the general case of high salt 
concentrations. We propose to take advantage of this 
similarity by applying insights gained from studies at 
high electrolyte concentrations to calculate the charge 
density at silica-water interfaces with low surface area 
and with no added salt. 

Whether the specific surface area is indeed small 
enough to warrant a "high salt" treatment, depends 
largely on the geometry of the considered experimen- 
tal setup and has to be checked on a case-by-case ba- 
sis. For some of the aforementioned interaction measure- 
ments 1 12 ItJ this approach is appropriate, in other cases 
it may provide a very rough upper limit for the surface 
charge. If, on the other hand, the counterions due to 
the charged surfaces give a non-negligible contribution 
to the overall ionic strength, they have to be considered 
explicitly, for instance within a cell model ■ Here, 

we concentrate on the former case of "high salt" and 
adopt the Basic Stern model |^] , which has been shown 
to accurately describe titration data |2^] obtained in the 
this regime for nominally nonporous, fully hydrated silica 
particles 

Within the Basic Stern model the charge of silica is 
regarded as localized entirely on the surface and arising 
from a concentration Fg i0 - of dissociated head groups 



SiO- 
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Under normal conditions, only a fraction of the total con- 
centration, 



r — r Si0 - + r S iOH, 



(3) 



of chargeable sites dissociate. The relevant mass action 
law for the deprotonation reaction, Eq. ([[]) , 
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is characterized by the logarithmic dissociation constant, 
pK, and accounts for the influence of the surface's elec- 
trostatic potential, ipo, through the surface activity of 
protons, 
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Here, [H+] b = l(T pH Mol/1 is the bulk activity of pro- 
tons, and = ksT denotes the thermal energy. 
The dissociation constant is an inherent property of the 
silicate- water interface and is estimated to be pK = 7.5 
on the basis of a surface complexation model |24[| . 

As counterions dissociate from the surface, they form 
a diffuse cloud of charge within the electrolyte. The 
Basic Stern model treats the counterions as being sepa- 
rated from the surface by a thin Stern layer across which 
the electrostatic potential drops linearly from is surface 
value, V'Oi to a value tpa called the diffuse layer potential 
|p8| , f22f . This potential drop is characterized by the Stern 
layer's phenomenological capacity, 
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This capacity, C, reflects the structure of the silicate- 
water interface and should vary little with changes in 
surface geometry or electrolyte concentration. Titra- 
tion data on colloidal silica |29) are consistent with C = 
2.9 F/m 2 j2|. 

Eqs. (0-|) can be solved for the diffuse layer potential 
as a function of the charge density on the interface: 
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This relation reflects the chemical nature of the interface 
and its charging process. 

Another functional dependence follows from the dis- 
tribution of mobile charges in the solution. If the latter 
is described by the Poisson-Boltzmann equation (PB), 
then the charge of an isolated, flat surface satisfies the 
Grahame equation 
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Here, e is the permittivity of the solution and the 
Debye screening length given by k 2 = 47r/3e 2 n/e, where 
n is the total concentration of small ions, all of which are 
assumed to be monovalent. The generalization of Eq. (||) 
to account for a curvature of radius a, 



27r/?e 



sinh 



tanh 



(9) 



is known to give the surface charge density to within 5% 
for na > 0.5 and any surface potential p0[|. 

Combining Eq. (R) with Eq. (||) or (S) yields self- 
consistent values for the surface charge density, <r, and 
the diffuse layer potential, ipd- These values character- 
ize the equilibrium of bound and mobile charges in the 
interfacial region, but are not necessarily accessible ex- 
perimentally, given the requirement of large surface areas 
for potentiometric titrations and the interpretive ambi- 
guities inherent to other techniques. 

Most measurements of interfacial interactions probe 
the electrostatic potential ip at distances for which eip < 
ksT. Under these conditions ip is described accurately 
by the linearized Poisson-Boltzmann equation, whose so- 
lution for a single flat surface has the form ip(x) = 
tpeft exp(— kx), where x is the distance from the surface. 
The effective surface potential ip e s in this experimentally 
accessible regime is related to the actual diffuse layer po- 
tential through pj| 



Peip e f[ = 4 tanh 
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Again, there is an approximate generalization for 
curved surfaces |3lJ]: 
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The associated effective charge density can be obtained 
from 
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which is just the linearization of Eq. (|^). This effective 
charge density characterizes essentially all of the recent 
measurements of electrostatic interactions between well- 
separated charged surfaces. 

The effective charge's relevance to experimental obser- 
vations is based in the popularity of the linear superposi- 
tion approximation for estimating the interaction energy, 
u{h), between two charged spheres of radii a\ and ai as a 
function of their surface-to-surface separation h. In this 
approximation |32||, 
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where <j\ and 02 should be understood to be effective sur- 
face charge densities obtained from Eq. (|ll|) rather than 
the bare charge densities from Eqs. (0) and @. Using 
the effective surface charge densities implicitly accounts 
for overexponential decay of the electrostatic potential 
near the surfaces that follows from the nonlinearity of the 
Poisson-Boltzmann equation. The interaction between a 
sphere and a planar wall is obtained by taking the limit 
of one infinite radius in Eq. (|l3|). 

Fig. shows computed values for the bare and effective 
charges of a planar silica surface and a 1 /im-diameter 
silica sphere for pH values between 7 and the lowest 
pH compatible with an ionic strength of 1, 5, and 10 
/1M0I/I. These arc reasonable values for deionized water 
under usual experimental conditions. In addition to using 
C = 2.9 F/m 2 and pK = 7.5, we have further assumed a 
total site density of V = 8 nm~ 2 , a commonly cited lit- 
erature value for nonporous, fully hydrated silica ||. Al- 
though r could vary widely depending on surface prepa- 
ration, the degree of protonation is determined mostly by 
the electrostatic interactions among the small fraction of 
charged surface sites, rather than the large number of 
neutral sites that T accounts for, and so our results are 
quite insensitive to this parameter. This robustness val- 
idates our assumption that details in the structure of 
nonporous surfaces do not matter in the present context. 
Indeed, the top graph of Fig. |l| also should represent the 
charging properties of a polished glass surface. Note how- 
ever that our arguments do not apply to some types of 
silica that are believed to be very porous and contain a 
much higher charge Q , the largest part of which seems 
located in the porous volume [ p3| rather than on the sur- 
face. 

Fig. [I] demonstrates that the effective charge densi- 
ties in deionized solutions do not depend as sensitively 
on pH as their "bare" counterparts, but that a signifi- 
cant variation with ionic strength persists. The top part 
of the figure indicates that the value of 2000 effective 
charges per /1m 2 (= —0.32 mC/m 2 ) assumed by Squires 
and Brenner M] for a glass plate in contact with deion- 
ized solution of = 0.275 /im (i.e. an ionic strength 
of 1.2 x 10 -6 M) is very reasonable. The confirmation of 
this previously very uncertain value provides vital sup- 
port for their recent electro-hydrodynamic explanation of 
apparent attractions between like-charged particles near 
a single glass wall. 

In a study of the equilibrium interaction between few 
1.58 /im silica sphere at the bottom of a large glass 
container filled with deionized water, we found an ionic 
strength between 8.5 x 10~ 7 and 1.1 x 10~ 6 M and a 
charge density between 550 and 830 e//im 2 from a fit of 
measured interaction energies to Eq. ([l3]) jlT]]. Compar- 
ison with Fig. [l] shows that these charge densities are a 
little below our expectation for isolated spheres, but have 
the right order of magnitude. The remaining difference 
can be explained by the spheres' close proximity to the 
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FIG. 1. The bare (full lines) and effective (dashed curves) 
charge densities of a planar glass wall and a 1 micron silica 
sphere, assuming a density r = 8 nm -2 of chargeable sites, 
a pK value of 7.5 for the silanol dissociation, and a Stern 
capacity of 2.9 F/m 2 . 



bottom wall of the glass container, as we describe in the 
next section. 



III. CHARGE REGULATION OF ANIONIC 
SURFACES 

Many experimental techniques to measure colloidal 
forces (AFM, TIRM, and the surface force apparatus, 
for instance) provide good resolution at very short dis- 
tances. Often the surface separations h of interest are 
comparable to or smaller than the screening length and 
much smaller than the radii of curvature, a\ and ai-, of 
either surface. In this situation, the surfaces may be re- 
garded as locally flat, and the Derjaguin approximation 



f(h) = 2n^^W(h). 
ai + a,2 



(14) 



accurately expresses the magnitude f(h) of the acting 
force in terms of the interaction energy per unit area 
W(h) for two parallel (thick) plates of the same separa- 
tion. This reduces the interaction problem to finding the 
energy W{h) or the dividing pressure U(h) — —dWjdh 
resulting from a one-dimensional distribution of mediat- 
ing small ions. 

On the other hand, the superposition of the noninter- 
acting surfaces' electrostatic potential is no longer war- 
ranted; nor can reliable results be expected from a so- 
lution of the linearized Poisson-Boltzmann equation, as 
it only captures the case of weak potentials (atypical 
for strongly interacting surfaces) or, with renormalized 
charges, the asymptotic behavior for large separations. 

Instead, we consider the nonlinear PB equation, which 
in one dimension and for an excess of monovalent elec- 
trolyte ions reads 

d 2 * 



dx- 
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sinh ^(x), 
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where ^ = (3eip is the dimensionless electrostatic poten- 
tial and x the coordinate normal to the surfaces. Apply- 
ing this mean field formalism to more general electrolytes 
shall not be discussed here, since neglected ion correla- 
tions and ion-specific interactions with the surfaces tend 
to complicate the case of polyvalent ions. 

From Eq. ( |l5| ) it is clear that ^(x) is a convex function 
for if? < 0. The charge density on either surface (1 or 2) 
is given, according to Gaufi' law, by 

Cl,2 : 



4:Trf3ez dx l 2 ' ^ ^ 

where the derivative is taken at the surface with respect 
to its outward normal. It follows that between two nega- 
tively charged surfaces, has a maximum ^ m , which 
for identical surfaces lies exactly at the midplane. Choos- 
ing generally the position of this maximum as the origin 
of our coordinate system (i.e. ^(O) = 'I'm), we can ex- 
press the solution of Eq. @ as HQ 

f>{x) = + 21ncd(u|m), (17) 

with 
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(19) 



m = exp (24'm) , 

where cd(u\m) is a Jacobian elliptic function of argument 
u and parameter to fjBf . The derivative is 

d^f? / i/A _ha\ sn(w|m) 



= U 3 / 4 



dx V J cn(u|m) dn(u|m) ' ^ 

where sn(u|m), cn(u|m), and dn = cn/cd are again Ja- 
cobian elliptic functions of the argument u and parame- 
ter to given above. Efficient numeric implementations of 
these functions are readily available from mathematical 
libraries |Bq|. 
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Equal Surfaces 

We shall measure the separation h between the sur- 
faces by the distance between the head ends of the dif- 
fuse layer, i.e., ^>(h/2) — (3eipd- Evaluating Eqs. ( [l6| - 
|20| ) at x — ±/i/2 and combining them with the chemical 
boundary condition, Eq. (Q), provides an expression for 
the actual surface charge density a. The second bound- 
ary condition, dty /dx = at x = 0, is already implied in 
the solution for *(af), Eq. fll7|). Note that for <j(<ip d ,h) 
defined by Eqs. ( |l6| ) through (20), the long distance limit 
lim/i—Kx, crtyd, h) is given by Eq. (|J), which we have pre- 
viously used for isolated surfaces and weakly interacting 
surfaces in the superposition approximation. In general, 
a is a function of ipd as well as of the surface separa- 
tion and cannot be expressed analytically. As before, the 
electrostatic definition of cr, Eqs. ([t6]) and (po|), depends 
parametrically on the Debye length, while the chemical 
definition, Eq. (H), depends only on pH and the surface 
chemical parameters, T, pK, and C. 

Technically, the combination of Eqs. and (|l6[j2C| ) 
results in a single transcendental equation for the mid- 
plane potential ^ m , whose numerical solution provides 
a very convenient alternative to numerical integration of 
the differential equation ( |l5[ ) with nonlinear boundary 
conditions. 

Solving for \l/ m has the further advantage of immedi- 
ately yielding the electrostatic force per unit area 



II = nfcsT (cosh * m - 1) 



(21) 



i.e. the excess osmotic pressure of small ions at the mid- 
plane where the electric field and the associated Maxwell 
stress are zero. 



Dissimilar Surfaces 

The procedure described before may be applied to 
negatively charged surfaces other than glass or silica as 
long as the chemically imposed charge-potential relation 
<j(ifid) is modified to account for the surface properties of 
the considered material. Carboxylated latex for instance 
can be described in the same framework as silica, with a 
pK value of 4.9 for the dissociation of the carboxyl sur- 
face groups (COOH ^ COO~ + H+) and a large Stern 
capacity C (any value C ^ 10 amounting to a negligible 
potential drop \ipo — ipd\ across the Stern layer) p7fl . Sul- 
fate latex, on the other hand, may be considered as hav- 
ing a constant charge density (cr(V><i) = — er = const.), 
because the strongly acidic sulfate groups are fully dis- 
sociated in all relevant solution conditions. Fig, p shows 
the predicted (and experimentally confirmed p4| 37|,|38|] ) 
charging behavior of the aforementioned materials. The 
site density of T = 0.25 nm -2 chosen for both the sul- 
fate and the carboxyl surface lies in the typical range for 
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FIG. 2. The charge density of surfaces with silanol, car- 
boxyl, or sulfate head groups as a function of pH. Parameters 
for the silica-like surface are as in Fig. |l|. For both the sulfate- 
and the carboxyl bearing surface we have assumed a density 
r = 0.25 nm -2 of sites, all of which are constantly charged 
in the sulfate case; further parameters of the carboxyl surface 
are a large (infinite) Stern capacity and a dissociation pK of 
4.9. 



commercially available latex spheres and has also been 
cited as the density of carboxyl groups on the membrane 
of blood cells |3 . 

A way to evaluate the interaction between two dissim- 
ilar surfaces starts by applying the described method for 
equal surfaces separately to both materials. For each 
of these symmetric systems, one obtains the midplane 
potential ^ m and thus via Eq. ( p7| ) the full potential 
function ^(x) associated with any given separation be- 
tween equal plates. Since ^>(x) is already fully deter- 
mined by the value of \l/ m = ^(O) and the requirement 
d^ /dx\ x= o = 0, solutions ^f(x) of the Poisson-Boltzmann 
equation for different systems are identical if they corre- 
spond to the same \l/ m (i.e. the same pressure), the only 
difference being the surface separation h for which they 
occur in the two systems. A solution ^(x) associated 
with a separation hi in one system and with separation 
hi in the second system clearly serves as a solution in 
a mixed system with a surface of type 1 at x = —h\/2 
and a surface of type 2 at x = +/12/2. Moreover, the 
separation h = (hi + h^jl at which this solution occurs 
in the mixed system is unique, because the pressure is a 
monotonic function of separation in the symmetric sys- 
tems and can thus be inverted to give the two separation 
functions /ii(^ r m ) and ^2(^111) • Our strategy therefore 
consists of computing \l/ m for all separations of interest 
in the symmetric systems 1 and 2, finding the separa- 
tions hi(^ m ), h2(^ m ) by inversion, and finally inverting 
their arithmetic mean ^(^m) = [/ii(^m) + ^2('J / m)]/2 to 



5 



obtain \& m (/i) and all the ensuing properties of interest 
in the mixed system. 

Some results of this type are shown in Fig. |3[ where 
we have plotted the charge density of a glass or silica 
surface and the electrostatic pressure as it interacts with 
either its own kind or with a surface of the carboxyl or 
the sulfate type. At the chosen ionic strength of 1 mM 
and pH 6, the charge of the silica surface is seen to de- 
viate significantly from its value in isolation (horizontal 
line and Fig. |^) up to separations of several screening 
lengths. Moreover, the nature of the second surface also 
has a profound effect not only on the strength of the 



interaction, but also on the charging state of the silica. 
While all anionic surfaces will reduce the effective charge 
on silica upon approach, the rate at which they do so 
strongly depends on the amount and variability of their 
own charge. Neither of these dependencies are usually 
considered in the discussion of interaction measurements. 
A recent attempt to determine these ionization proper- 
ties of silica experimentally with atomic force microscopy 
p9[ has been limited to symmetric surfaces, and relies on 
model assumptions both for the charge regulation and for 
the strong van der Waals forces at short surface separa- 
tions. 
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FIG. 3. The charge density of the silica surface and the 
force per unit area it experiences when interacting with any 
of the three presented types of surfaces at pH 6 and an ionic 
strength of 1 mM = 9.6 nm). Surfaces parameters as in 
the previous figures. 



IV. CONCLUSIONS 

We have seen that the apparent charge on glass and 
silica surfaces of low specific area in pure water can be 
understood in terms of a simple model that was originally 
developed and tested for high electrolyte concentrations. 
Model predictions for effective charge densities compare 
favorably with interaction experiments on highly diluted 
silica spheres in deionized water fT^ ]. They also sup- 
port a new kinematic explanation of spurious long range 
attractions between like-charged particles near a single 
glass wall Jl6j . 

The regulated charge of silica and glass surfaces near 
contact with a second anionic surface, as well as the 
strength of the interaction, has been calculated from an 
exact solution of the nonlinear Poisson-Boltzmann equa- 
tion. In commonly encountered solution conditions, the 
charge regulation of silica was found to be effective at 
separations well beyond a Debye length. It also proves 
very sensitive to the chemical nature of the opposing sur- 
face. Although the additional presence of van der Waals 
forces makes a quantitative measurement of these effects 
difficult, they should certainly be accounted for in the 
interpretation of interaction experiments. 

This work was supported by the National Science 
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by the Deutsche Forschungsgemeinschaft. 
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